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ABSTRACT 

Formation of interstellar clouds as a consequence of thermal instability is studied using two- 
dimensional two-fluid magnetohydrodynamic simulations. We consider the situation of converging, 
supersonic flows of warm neutral medium in the interstellar medium that generate a shocked slab 
of thermally unstable gas in which clouds form. We found, as speculated in paper I, that in the 
shocked slab magnetic pressure dominates thermal pressure and the thermal instability grows in the 
isochorically cooling, thermally unstable slab that leads formation of HI clouds whose number density 
is typically n < 100 cm -3 , even if the angle between magnetic field and converging flows is small. 
We also found that even if there is a large dispersion of magnetic field, evolution of the shocked slab 
is essentially determined by the angle between the mean magnetic field and converging flows. Thus, 
the direct formation of molecular clouds by piling up warm neutral medium does not seem a typical 
molecular cloud formation process, unless the direction of supersonic converging flows is biased to the 
orientation of mean magnetic field by some mechanism. However, when the angle is small, the HI 
shell generated as a result of converging flows is massive and possibly evolves into molecular clouds, 
provided gas in the massive HI shell is piled up again along the magnetic field line. We expect that an- 
other subsequent shock wave can pile up again the gas of the massive shell and produce a larger cloud. 
We thus emphasize the importance of multiple episodes of converging flows, as a typical formation 
process of molecular clouds. 

Subject headings: magnetohydrodynamics — ISM: general — method: numerical — instabilities 



1. INTRODUCTION 

Molecular clouds are the sites of all the present-day 
star formation. However, our understanding of the 
physical conditions of molecular clouds is very limited. 
One promising method to reveal the physical conditions 
of molecular clouds is studying formation of molecu- 
lar clouds. As a formation mechanism of molecular 
clouds, converging, supersonic HI flows that trigger ther- 
mal instability to form cold gas clouds have been studied 
(Ballesteros-Paredes et al. 1999; Hartmann et al. 2001; 
Inoue & Inutsuka 2008; Hennebelle et al. 2008; Baner- 
jee et al. 2008; Heitsch et al. 2005, 2006, 2008, 2009; 
Vazquez- Semadeni et al. 2006, 2007, see also Koyama & 
Inutsuka 2000, 2002). 

Recently, using magnetohydrodynamic (MHD) simu- 
lations, Hennebelle et al. (2008), Banerjee et al. (2008) 
and Heitsch et al. (2009) performed the formation 
of molecular clouds essentially considering a situation 
where the orientation of the converging flows and mean 
magnetic field are parallel. However, in Inoue & Inut- 
suka (2008; henceforth paper I), we showed in a two- 
dimensional two-fluid MHD simulation that, if the con- 
verging flows and magnetic field are orthogonal, only the 
formation of HI clouds (n < 100 cm -3 ) is possible instead 
of the formation of dense clouds (precursor of molecular 
clouds; n > 100 cm -3 ) due to the effect of magnetic pres- 
sure. Using three-dimensional MHD simulation, Heitsch 
et al. (2009) also obtained results similar to ours even in 
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the case with weak magnetic field (B = 0.5 fiG). 

In the converging flow scenario, to trigger cooling col- 
lapse or thermal instability, supersonic flows are neces- 
sary (Hennebelle & Perault 1999). Since magnetic pres- 
sure in the diffuse ISM is comparable to thermal pres- 
sure, it is natural to assume that the direction of super- 
sonic flows and mean magnetic field has no strong corre- 
lation. Therefore it is important to investigate whether 
the formation of molecular clouds is possible when the 
converging flows and mean magnetic field have a small 
angle. Using one-dimensional MHD simulation, Hen- 
nebelle & Perault (2000) showed that if the initial an- 
gle between the magnetic field and the fluid velocity is 
larger than 20-30 degrees with a few microGauss mag- 
netic field strength, magnetic pressure prevents growth 
of thermal instability. In multi-dimension, thermal in- 
stability can always grow along the magnetic field line 
(paper I). But the effect of magnetic pressure drastically 
change the physical condition of clouds formed by ther- 
mal instability, such as density and temperature (paper 
I). In this paper, using two-dimensional MHD simula- 
tions, we extend paper I and examine the cases where the 
converging flows and mean magnetic field have angles as 
Hennebelle & Perault (2000). This paper is organized as 
follows. In §2, we provide numerical setup of simulations. 
The results of simulations are shown in §3. Finally, in §4, 
we summarize our results and discuss their implications. 

2. NUMERICAL SETUP 

As in paper I, we consider a weakly ionized interstel- 
lar medium using two-fluid approximation incorporating 
with the effects of coupling between neutral and ionized 
gases ( friction force/heating, ionization, and recombina- 
tion), radiative cooling/heating, and thermal conduction. 
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Details of basic equations (ideal MHD and HD equations 
with the source terms of two- fluid coupling, heating, cool- 
ing, and thermal conduction) and numerical schemes to 
solve them were developed in paper I. In Fig. HJ the ther- 
mal equilibrium state of the cooling function we used is 
shown as a thick solid line. We also overlay the dotted 
regions in which a gas is thermally unstable (see, Bal- 
bus 1995; Field 1965; Schwarz et al. 1972; Koyama & 
Inutsuka 2002 for linear stability analyses). 

We use a two-dimensional, rectangular, numerical do- 
main whose side lengths are L x = 30 pc and L y = 10 
pc that is divided by N x x N y uniform square cells 
(N x = 3N y ). As an initial condition, we set converging 
flows of warm neutral medium (WNM) that collide head- 
on at the center of x-axis. Density and thermal pressure 
of the WNM are n = 0.67 cm -3 and p th /k B = 4,000 
K cm -3 , which are in thermal equilibrium and typical 
values of WNM in the ISM. Our initial condition for the 
ionized gas is ionization equilibrium to the neutral gas. 
Interested readers should refer to Paper I for details. The 
velocity field of the converging flows is given such that 



v x (x,y)=v (x) + 5v x (x,y), 
v y {x,y) = 5v y {x,y), 



(1) 
(2) 



where vq(x) = v cnv > (x < L x /2) and = — v cnv < 
(x > L x /2) is the converging component of the velocity, 
and Sv Xi y is the fluctuation component. The fluctuation 
component is given by 
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where 0i 5 2 are random phases and function A is chosen so 
that fluctuations show the Kolmogorov energy spectrum 
E(k) oc /c -5 / 3 and velocity dispersion of the fluctuations 
is set at 10% of the sound speed of the WNM. Such a 
spectrum and velocity dispersion are naturally expected 
in the ISM, since WNM is turbulent with velocity dis- 
persion being transonic at a scale of a few hundreds pc 
(Dib et al. 2006), which leads velocity dispersion roughly 
10 % of the sound speed at scale of a ten pc. A uniform 
magnetic field with an angle of to the x-axis is ini- 
tially imposed. According to the observations using ro- 
taion/dispersin measures (Rand & Kulkarni 1989; Beck 
2000), the magnetic field strength in the Galactic plane is 
about 1-5 fiG (approximately equipartition to the ther- 
mal energy). Thus, we use the fixed strength of mean 
magnetic field Bo = 3.0 fiG throughout of this paper 
that gives plasma (3 = 1.5 and meets the above observa- 
tions. 

We also perform simulations with an initially tangled 
magnetic field, instead of the velocity fluctuations. In 
this model, the magnetic field is given by 



B x (x,y)=B : 
B y (x,y)=B ] 



mean 
mean 



COS(0) + ^tangle,^, 2/), (5) 

sin(6>) + £ tang ie, y (x, y), (6) 



TABLE 1 
Model Parameters 



Name v C nv [km s 



[deg.] 



B n 



(-^tangle )r 



la 20.0 


15.0 


Bo 


= 3.0 /j,G 


N/A 




lb 20.0 


40.0 




Bo 


N/A 




2a 10.0 


15.0 






NT / A 
IN / A 




2b 10.0 


40.0 




Bo 


N/A 




3a 20.0 


15.0 




Bo 


Bo 




3b 20.0 


15.0 




Bo 


2.0 S 
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-^tangle, x (^j V) ~ 


d y A z (x,y), 






(7) 


-^tangle, x V) ~ 


&x A z (x, 


y), 






(8) 


10 

A z (x,y)= Y C(k) cos 


"2tt 


(k x x + k y y) 


+ 0(fc) 


(9) 


k x iky - 


=0 











where the tangled components of the magnetic field are 



We use random phase as 4>{k). The spectrum of the 
magnetic field in the ISM is highly unknown. Since the 
spectrum of the magnetic field fluctuation is so unknown, 
we use C(k) oc k~ 4 following Heitsch et al. (2009). Thus, 
this model of magnetic field is characterized by the mean 
magnetic field strength 5 mean , the angle and a root- 
mean-square of the fluctuating field (5 ta ngie)rms- 

We impose periodic boundary conditions at y = and 
L y surfaces. The fluctuating converging WNM flows are 
set at x = and L x surfaces. The boundary condition 
of the fluctuation component is given by substituting x 
for — t (at x = 0) and for L x — vo t (at x = L x ) in eqs. 
®, Q, and ©. 

We examine various situations by changing the free 
parameters of the models, which are summarized in Table 
1. In this paper, we concentrate on the cases where not 
equal to and 90 deg., more specifically 9 = 15 and 40 
deg., since the case of 6 = is examined in Hennebelle 
et al. (2008), Banerjee et al. (2008), and Heitsch et al. 
(2009) and the case of = 90 is examined in paper I and 
Heitsch et al. (2009). According to the one-dimensional 
study by Hennebelle & Perault (2000), the evolution of 
the shocked medium is affected greatly by the magnetic 
pressure, when the angle is larger than 20-30 deg. with 
a few microGauss initial field strength. Thus, we choose 
the cases 6 = 15 and 40 deg. that are respectively smaller 
than and larger than the angle 20-30 deg. 

Our initial conditions may correspond to the propaga- 
tion of shock waves in the ISM, e.g., super-shell due to 
multiple supernovae and galactic spiral shock. In general, 
spiral shocks are highly oblique. However the oblique 
shock can be transformed to a simple one-dimensional 
shock by the appropriate Galilean transformation. Thus, 
our analyses with various angles of magnetic field lines 
with respect to the shock normal can also be applied to 
spiral shocks. In the case of the super shell, one of the 
ram pressures of converging flows would be substituted 
by the thermal pressure of the coronal gas, and in the 
case of the spiral shock that would be substituted by 
enhanced thermal pressure due to stellar potential. In 
super-shell, after several million years from expansion at 
which thermal instability is effective in the shell, prop- 
agation speed of shock wave is a few times ten km/s 
depending on the number of SNe (see, e.g., Tomisaka 
1998). Thus, we investigate the cases ^cnv — 
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number density : n [ cm 3 ] 

Fig. 1. — Thermal equilibrium state of our cooling function (thick 
line). The dashed lines represent isothermal lines of T = 10 4 , 10 3 
and 10 2 K. In the dotted region, gas is thermally unstable. The 
typical gas evolution of Model la and lb is drawn schematically as 
arrows. 

km/s. 

As shown in the following section, the results of most 
models show formation of clouds whose density is n < 
100 cm -3 . In such cases the Field length that char- 
acterizes the smallest scale of the system is approxi- 
mately > 0.01 pc (Inoue et al. 2006). Thus, in or- 
der to resolve the Field length, we use the resolution 
N x x N y = 3072 x 1024. The data of the computations 
are saved each 0.2 Myr of the simulations that are used 
in the analyses given below. 

3. RESULTS 
3.1. Results of Model la and lb 

We show density structures and magnetic field lines 
of Model la and lb in Fig. [2] (a) and (b), re- 
spectively. The depicted times are chosen when 
the thickness of shocked slab becomes larger than 
20 pc. Here, we define the thickness of the slab 
such that / s i a b = (£ s , r i g ht (y)) y {x 8 ,\eft(y)) y , where 
(•••}y means average over ^/-direction, £ S) i e ft(y) = 
min[x | B y (x,y) 2 /8ir > 3 x 10~ 13 g cm 2 s~ 2 ], and 
fright (2/) = max[x | B y (x, 2/) 2 /8tt > 3 x 10 13 g cm 2 s~ 2 ]. 
In the two models, structures of shocked slabs are simi- 
lar. The evolution of the shocked WNM can be under- 
stood as follows (see also paper I): At first, the WNM is 
compressed by the shock by which the WNM kicked out 
from the thermal equilibrium states and into states dom- 
inated by cooling. Then the slab of shocked cooling gas 
contracts nearly isobarically in order to oppose the ram 
pressure of the converging flows. This contraction leads 
amplification of magnetic field by compressing transverse 
components of the magnetic field. Once the magnetic 
pressure of the amplified magnetic field becomes compa- 
rable to the ram pressure, the isobaric contraction stops 
and the gas begins to cool nearly isochorically. The iso- 
baric phase of the evolution does not last long, because 
the transverse component of the magnetic field is already 
amplified by the shock compression and the contraction 
by only a factor of a few can strengthen the magnetic 
field enough to oppose the ram pressure. During the iso- 
choric evolution phase, the gas becomes thermally un- 
stable and some portion of the gas, whose density is en- 
hanced by the compressions due to the turbulence given 



Time = 6.60 Myr 




5 10 15 20 25 30 

Time = 6.40 Myr 




5 10 15 20 25 30 



x[pc] 

Fig. 2. — Structures of number densities and magnetic field lines. 
Field lines are depicted so that distance from the neighbor field 
line lyn is inversely proportional to the field strength (B = 10 /iG 
[Z m /1.5pc] -1 ). Panels (a)-(f) represent the results of Model la, 
lb, 2a, 2b, 3a and 3c, respectively. 

initially, grows toward a stable dense equilibrium, i.e., 
cold neutral medium (CNM). Since thermal instability 
evolves essentially along the magnetic field, filaments of 
CNM are formed. On the other hand, the remaining gas 
continues to be cooled isochorically, until it reaches un- 
stable equilibrium and then grows toward a stable warm 
phase (WNM). The sequence of the typical evolutions is 
drawn schematically in Fig. [TJ 

In Fig. [3] (a) and (b), we plot the probability dis- 
tribution function (PDF) of the gas that confirms the 
evolution discussed above. Note that since the timescale 
of thermal instability toward the WNM phase is approx- 
imately 5 Myr (Inoue et al. 2006; 2007), which is longer 
than the timescale of the instability toward the CNM 
(< 1 Myr), the diffuse gas component in the shocked 
slab is still growing toward the WNM phase. In Fig. [4] 
(a) and (b), we plot cross section profiles of the thermal 
(pth)? magnetic (By/Sn), total (p t h -\- By/Sir), and ram 
pressure (pv%) at y = 5.0 pc. We can confirm that the 
shocked slab is supported by the magnetic pressure. Note 
that From the instant PDF in Fig. [3] and cross section 
in Fig. [U we can also know the information about time 
evolution of shocked gas. The gas in the downstream at 
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Fig. 3. — PDFs in n—P plane. Panels (a)-(f) represent the 
results of Model la, lb, 2a, 2b, 3a and 3c, respectively. Thermal 
equilibrium state (solid line), initial point of simulation (cross), 
and isotherms (dashed lines) are plotted. 



the distance r from the shock front roughly shows the 
state of the gas at the time r/vj since shocked, where 
Vd is the downstream velocity at the shock rest frame. 
From this, we can also see that the magnetic pressure 
promptly dominates the thermal pressure inside the slab 
that corresponds to the briefness of the isochoric con- 
traction phase. 

Densities of the CNMs are typically < 100 cm -3 , and 
temperatures are ~ 50 K corresponding to HI clouds. 
High density clumps of n > 100 cm -3 are formed only at 
an early stage around the center of the numerical domain 
(x ~ 15 pc), since high density CNMs can be generated 
until the thermal pressure of the slab begins to decrease, 
i.e., until the magnetic pressure dominates the thermal 
pressure. To see this, we plot the average column density 
evolutions of the high density (n > 200 cm -3 ) and low 
density (200 cm _3 > n > 20 cm -3 ) CNMs along the pr- 
axis in Fig. [5] (a). 

Note that only the magnetic field amplification a few 
times the initial strength makes the opposition of the 
magnetic pressure and the ram pressure possible. Thus, 
the magnetic field strength in the slab is not so high in the 
simulations (B ~ 10 fiG) that will decrease to average 
field strength in the ISM once the converging flows stop. 

3.2. Results of Model 2a and 2b 

We show density structures and magnetic field lines of 
Model 2a and 2b in Fig. [2] (c) and (d), respectively. As 
in Model la and lb, the depicted times are chosen when 
the thickness of the shocked slab becomes larger than 20 
pc. Since the velocity of the converging flows is smaller, 
the shock waves propagate faster than in Model la and 
lb. 

The two models are apparently similar. However, de- 
tails of structures and their evolutions are different. In 
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Fig. 4. — Cross section profiles of thermal (red), magnetic 
(green), total (blue) and ram (magenta) pressures at y = 5.0 pc. 
Panels (a)-(f) represent the results of Model la, lb, 2a, 2b, 3a and 
3c, respectively. 

the case of Model 2a four shock waves are induced, i.e., 
fast and slow shocks propagate to ± x directions, while 
slow shocks are not clearly seen in Model 2b like Model 
la and lb. This difference is seen more clearly in ear- 
lier stages. In Fig. [H we show density structures and 
magnetic field lines of Model 2a and 2b when the thick- 
ness of the shocked slabs is ~ 10 pc. In general, fast 
shocks amplify the magnetic field by compressing field 
lines, while slow shocks lower field strength by stretching 
the field lines, i.e., in the case of fast shock, the angle be- 
tween magnetic field and shock normal becomes larger in 
downstream compared to the upstream, and in the case 
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Fig. 5. — Average column density evolutions of the high density 
(n > 200 cm -3 ) and low density (200 cm -3 > n > 20 cm -3 ) 
CNMs along the x-axis. Panels (a)-(c) represent the results of 
Model la, lb, 2a, 2b, 3a, and 3b, respectively. Thick (thin) lines 
show the column density evolutions of high (low) density CNMs. 



of slow shock, the angle becomes smaller. By using these 
basic characteristics of MHD shocks, we can identify in 
the top panel of Fig. 6 (Model 2a) that there are fast 
shocks at x c± 10 and 20 pc and there are slow shocks at 
x c± 12.5 and 17.5 pc. On the other hand, in the bot- 
tom panel of Fig. 6 (Model 2b), there are only two fast 
shocks. 

The structures of the pressures in the two models also 
differ. Fig. [H (c) and (d) shows the cross section profiles 
of the thermal, magnetic, total (thermal plus magnetic), 
and ram pressures at y = 5.0 pc. In Model 2a, ow- 
ing to the slow shocks magnetic pressure does not dom- 
inate thermal pressure in the shocked slab. In this case, 
the condensation driven by cooling can proceed with- 
out decreasing thermal pressure greatly (nearly isobaric 
contraction) that leads formation of dense CNM whose 
number density is larger than 200 cm -3 . This evolution 
is similar to the evolution of shocked gas in unmagne- 
tized media (see, paper I). In the case of Model 2b, on 
the other hand, the magnetic and thermal pressures are 
comparable (see Fig. IH[d]). 

The magnetic pressures, in particular those in Model 
2b, are still increasing to dominate the thermal pres- 
sure. Since the shock propagation speed is fast, the 
shock waves reach computational boundaries within a 
few Myr that makes the timescale short within which we 
can examine the evolution of the shocked slab. Thus, we 
have not confirmed the clear domination of the magnetic 
pressure in the slab. However, in Model 2b, we expect 
from the absence of the slow shocks that the shocked slab 
would evolve into a two-phase medium composed of low- 
density CNM (n < 200 cm" 3 ) and WNM as in Model 
la and lb. In order to confirm the domination of mag- 




FlG. 6. — Structures of number densities and magnetic field lines 
when the thickness of the shocked slabs is ~ 10 pc. Field lines are 
depicted six times thicker than Fig. [2] Top and bottom panels 
show the results of Model 2a and 2b, respectively. 

netic pressure, simulation with a much larger numerical 
domain is necessary. 

We can also see the difference between the two models 
in Fig. 03(b). In Model 2a, the growth of column density 
of high density CNM does not show saturation unlike 
Model la and lb, while Model 2b evolves like Model 
la and lb, although the timescale we have examined is 
short. Note that the PDFs in n-P plane of Model 2a 
and 2b plotted respectively in Fig. [3] (c) and (d) does not 
clearly differ between the two models. This is because the 
slab of Model 2b is still in the isobaric contraction phase. 
The difference would be clear once the magnetic pressure 
dominates the thermal pressure in Model 2b that leads 
the evolution of the isochoric cooling in n-P plane. 

What controls bifurcation of the number of induced 
shocks? The condition of the bifurcation can be derived 
from the MHD Rankine-Hugoniot relations: 



s(po ~ Pi)=Po v x,o - Pi v Xi i, 



(10) 



S(P0 %0 - Pi %l) 



P0 V Xy 0V y ,0 - 

Pi v Xi i v V) i 



4ir 
B x B V) i 

47T 



,(H) 

(12) 



s (By,o — By^i) = (v Xj o B y ^ — v Vj o B x ) 

— ( V X,1 By^ — Vy^ B x ), 

where the subscripts and 1 respectively denote the up- 
stream and downstream value, s is the shock propagation 
speed, and we have used the fact that B x = constant 
across the shock wave owing to V • B = 0. The equations 
(l)-(3) represent the conservation laws of mass, tangen- 
tial component of momentum, and tangential component 
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Fig. 7. — Critical velocity as a function of with the typical 
parameters B = 3.0 /iG, no = 0.67 cm -3 , and r = 4. Conditions 
of Model la, lb, 2a, and 2b are plotted as crosses. 

of electric field, respectively. 

Immediately after the collision of the converging flows, 
v Xj i and ivo should be zero. By using these conditions 
and eqs. ([lQ |) -([T2 ]h eliminating s and v Vi i, we obtain 



(p 1 -p )B% Pov x , 



Po Pi v x ,o 



pi - po (By }0 - B y ,i) 



(13) 



In the case of the fast shock, if B y ^ > (< 0), B y $— B y ^ 
should be less (greater) than zero. Therefore, collision of 
the converging flow induces only the fast shock, when 



v x ,o > 



1 B cosO 



V 47r Po 



(14) 



where we have introduced the compression ratio of the 
density r = pi/po, and is the angle between converging 
flows and upstream magnetic field. If the velocity of 
the converging flows is smaller than the critical velocity 
^cnv < v CT , additional slow shocks are necessary in order 
to obtain physical down stream state. In Fig. [7J we show 
the critical velocity as a function of by using the typical 
parameters B = 3.0 /iG, no = 0.67 cm -3 and r = 4. We 
also plot the conditions corresponding to Model la, lb, 
2a, and 2b. 

In order to form dense CNM (precursor of molecular 
clouds) whose number density is larger than 200 cm -3 , 
the shocked gas should condense without decreasing ther- 
mal pressure, i.e., the magnetic pressure should not dom- 
inate in the slab. This requires the slow shocks that lower 
magnetic field strength in the slab. Thus, if the velocity 
of converging flows is smaller than the critical velocity a 
dense CNM sheet is formed as Model 2a. On the other 
hand, if the velocity of converging flows is larger than the 
critical velocity, the magnetic pressure dominate in the 
slab that leads isochoric cooling and formation of low- 
density HI clouds (n < 100 cm -3 ) as Model la, lb, and 
as expected in Model 2b. 

Therefore, only the case where v cnv < v cr (Model 2a) 
can generate molecular clouds directly from WNM. How- 
ever, the column density of dense CNM sheet formed in 
Model 2a (the precursor of molecular cloud) is at least an 
order of magnitude smaller than that of molecular clouds 
(> 10 21 cm -2 ). Thus, in order to form molecular clouds, 
the converging flows have to last a long time. Note that, 
even in such case, the shocks continue to propagate out- 



wards making direct formation of molecular clouds dif- 
ficult as follows: The shocked gas that accretes to the 
dense CNM sheet is provided with sub-slow speed far 
from the dense CNM sheet. In the ISM, such a slow flow 
is easily affected by, e.g., SNe and stellar winds by which 
flow direction would change. As a result, the accretion to 
the dense CNM is stopped leading to failure of molecular 
cloud formation. 

3.3. Results of Model 3a and 3b 

In Model 3a and 3b, instead of the velocity fluctua- 
tions, we impose the highly tangled magnetic field. Fig. 
[2] (e) and (f) respectively show density structures and 
magnetic field lines of the results of Model 3a and 3b at 
which the thickness of the shocked slab becomes larger 
than 20 pc 3 . We see that the results are close to that of 
Model la. The PDFs in Fig. [3] (e) and (f), the cross- 
section profiles in Fig. [H (e) and (f), and the average 
column density evolutions in Fig. [5] (c) also show be- 
havior similar to those of Model la and lb. Since the 
condition in Model 3b is more turbulent than Model 3a 
due to the larger tangled component of magnetic field, 
the distribution of the gas in n-P plane in Model 3b is 
more dispersed than Model 3a. 

From these results, we conclude that evolution of the 
shocked slab is essentially determined by the angle be- 
tween the mean magnetic field and converging flows, 
even if the mean magnetic field is weaker than the tan- 
gled component. The reason why the mean magnetic 
field determines the evolution of the shocked slab is that 
the tangled components of the magnetic field behave as 
Alfven waves which eventually decay due to the cascade 
of MHD turbulence within Alfven crossing time (Goldre- 
ich & Sridar 1995; Cho & Lazarian 2005). Therefore, if 
Vcnv > v cr , even when the tangled components dominate 
the mean field, the fate of the shocked slab is determined 
by the mean magnetic field. 

One may wonder the different conclusion given in 
Heitsch et al. (2009) in which tangled magnetic field 
component affect the post shock structure greatly. In 
Heitsch et al (2009), they compared the models in which 
the mean magnetic field is parallel to the flow direction 
except their Model Y05. In these cases, the shocked gas 
evolves in different ways depending on the local orien- 
tation of magnetic field as follows: If the local mag- 
netic field is parallel to the converging flows, the mag- 
netic field is not enhanced that would leads formation of 
dense CNM. On the other hand, if local magnetic field is 
perpendicular to the flow, the enhanced magnetic pres- 
sure would prevents formation of dense CNM that leads 
void. Thus, the tangled component of magnetic field can 
change the morphology of the clouds, if we consider the 
case of vanishing transverse mean magnetic field. In our 
calculation, we consider the models in which the mean 
magnetic field and the converging flows have finite angles. 
In these cases, if v cnv > v cri the magnetic pressure even- 
tually dominate the thermal pressure in the post shock 
media that prevents direct formation of dense CNM and 
leaves two-phase media composed of low-density CNM 



3 In Model 3b, 
min(max)[ x | B y (x : y) 2 /8tt 



s,left(right) (y) 
21 



we define 

> 3 x 10 -12 g cm 2 s ~" j, since 

the initial transverse component of magnetic field in Model 3b is 
larger than the other models. 
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and WNM. 

4. SUMMARY AND DISCUSSION 

4.1. Summary 

In this paper, using two-dimensional two-fluid MHD 
simulations, we studied the formation of interstellar 
clouds as a consequence of thermal instability in the 
shocked slab generated by converging WNM flows. Our 
findings are as follows: 

• The evolution of the shocked slab is controlled by 
the parameters v cnv (the velocity of the converg- 
ing flows) and (the angle between the converging 
flows and the mean magnetic field). 

• If the velocity of the converging flows is larger than 
the critical velocity v CY derived in eq. (jHj) and 
is at least larger than 15 deg., HI clouds embedded 
in WNM are formed in the slab, as we speculated 
in paper I. In general, in the case with a moderate 
magnetic field strength (B = 3 /iG), the critical ve- 
locity v cr is comparable to the sound speed of the 
WNM ~ 10 km s _1 (see, Fig. g]) depending on 0. 
It is known that in order to trigger cooling collapse 
or thermal instability, supersonic flows are neces- 
sary (Hennebelle & Perault 1999). Therefore, this 
would be the most probable fate of the converging 
flows in the ISM. 

• If the velocity of the converging flows is smaller 
than the critical velocity or the angle is nearly 
zero, dense clouds that possibly evolve into molec- 
ular clouds can be formed in the slab (see, Hen- 
nebelle et al. 2008; Banerjee et al. 2008; Heitsch 
et al. 2009 in the case of 6 = 0). 

• The entanglement of the magnetic field is not an 
important control parameter for the evolution of 
the shocked slab. The important parameter is the 
angle between mean magnetic field and the con- 
verging flows. Even if strength of the mean mag- 
netic field is weaker than the strength of the entan- 
gled component (dispersion), the mean magnetic 
field determines evolution of the shocked slab. 

In summary, the direct formation of molecular clouds 
from WNM by converging flows is possible, only if the 
angle between the mean magnetic field and the converg- 
ing flows is close to zero or the velocity of converging 
flows is less than the critical velocity given by eq. [Ml 

Note that our results are obtained using two- 
dimensional simulations that may change the configu- 
ration of clouds in three dimensions due to the modes 
of fluid motion that are not appear in planar case, e.g., 
the interchange mode. Thus, 3D simulation should be 
performed in order to confirm our results. However, the 
model Y05 (0 = 90 deg. and B = 0.5 /iG) of Heitsch 
et al. (2009) shows that, even when the initial magnetic 
field strength is an order of magnitude weaker than the 
equipartition, the magnetic pressure dominate thermal 
pressure in the shocked slab. This suggests that there is 
no mechanism that reduces the enhanced magnetic pres- 
sure even in 3D. Thus, we can expect that the enhance- 
ment of magnetic pressure makes the direct formation of 
molecular cloud difficult even in 3D. 



In this paper, we also neglect the effect of self-gravity. 
As discussed in Koyama & Inutsuka (2002), the scale of 
the clumps of CNM generated by thermal instability is 
much smaller than the Jeans length. Thus, we can ne- 
glect the effect of self-gravity on each individual CNM 
clumps. Furthermore, the mean densities in the shocked 
slabs are smaller than 10 cm -3 (see Fig. 8 below) indi- 
cating the Jeans length of shocked slab larger than 50 pc 
and the free-fall time larger than 50 Myr, even if we un- 
derestimate the mean temperature of the shocked slabs 
as 100 K. Thus, we can safely neglect the effect of self- 
gravity that leads contraction of shocked slab to form 
molecular cloud. Note that this is not the case when 
theta = deg, since the mean density of the slab would 
be much larger. 

4.2. Discussion 

We have shown that the direct formation of molecular 
clouds from WNM seems difficult unless some mecha- 
nisms creates biased supersonic flows of WNM that is 
parallel to the mean magnetic field orientation. Our re- 
sults would be also applied to the situations of propaga- 
tion of single shock wave, e.g., super-shell due to multiple 
super novae (SNe) and spiral shock wave. In the cases of 
the super shell, one of the ram pressures of converging 
flows that bound the shell would be substituted by the 
thermal pressure of the hot gas, and in the case of the 
spiral shock that would be substituted by enhanced ther- 
mal pressure that balances the stellar potential. 

However, if we consider some mechanisms that further 
pile up the gas in the two-phase shell generated by the 
converging flows, the formation of molecular clouds is 
possible. For example, in the case of Model la, the mean 
density of the gas in the two-phase shell is ~ 10 cm -3 
that is more than ten times higher than WNM. Thus, 
molecular cloud with Ay ~ 1 mag. would be formed if 
the two-phase gas in the shell is piled up again only over 
tens of parsec along the magnetic field line. The ISM is 
known to be frequently (once in a million years) swept 
up by SN shocks, which propagate a few hundred parsecs 
until decay. Thus, it can induce formation of molecular 
clouds by piling up the massive two-phase shell. As in 
the case of the converging flows studied in this paper, 
SN shock can pile up gas only when it propagates along 
the magnetic field line in the two-phase shell. Despite 
this, it is promising because of the high event rate of 
SN. We note the difficulty in propagation of a single SN 
shock wave piling up enough WNM to generate molecular 
clouds even if it propagates along magnetic field line be- 
cause of its short lifetime. The single SN shock can gen- 
erate molecular clouds only when the density of preshock 
gas is enhanced by something just as we consider here. 

What is the angle 6 below which the converging flows 
form the massive shell that can evolve into molecular 
clouds? If we assume an SN shock wave that pile up 
again the two-phase shell, the mean density in the two- 
phase shell should be larger than 

<n> > 6.3 cm" 3 ( -^-) ( , (15) 

where Z$n is the propagation length of SN shock until 
decay and Ay is the visual extinction of the formed 
molecular cloud. In §3.1, we have shown that the 
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Fig. 8. — Solid line shows the theoretical post shock mean 
density derived in equation (16). The parameters po and v x ^ are 
the same to Model la and lb and we use B y ^{6) = Bcos(#), where 
B = 3 fiG. The post shock values m and B Vi i are calculated by 
solving MHD shock jump conditions using the preshock conditions 
of Model la and lb. The mean densities of the results of Model la 
and lb are also plotted as crosses, which are evaluated at the times 
depicted in Fig. [2] The dotted line represents the critical mean 
density [eq. (15)]. If the mean density is larger than the critical 
value, a single SN shock can generate molecular cloud. 

mean density enhancement due to the isobaric contrac- 
tion of the shocked gas is stopped when the ram pres- 
sure of the converging flows and the magnetic pressure 
in the shell are balanced. This indicates that the iso- 
baric contraction enhances the post shock gas by a fac- 
tor of (pov^o)/(By^/87r), where subscript stands for 
preshock state (WNM) and subscript 1 stands for post- 
shock state (immediately behind the shock front). Note 
that B Vj i depends on (the angle between upstream 
magnetic field and shock normal). Thus, the mean den- 
sity of the shocked slab (two-phase shell) can be esti- 
mated by 



(n) ~ rii 



B 2 yA (0)/S^ 



(16) 



where subscript stands for preshock state (WNM), sub- 
script 1 stands for postshock state (immediately behind 
the shock wave), the x- and ?/-axes respectively represent 
directions normal to and tangential to the shock front, 
and is the angle between the preshock magnetic field 
and flow. In Fig. El we plot the mean density (n) (eq. 



[16]) against the angle between the converging flows and 
mean magnetic field 6 by using preshock conditions of 
Model la and lb. The postshock condition is calculated 
by solving the MHD Rankine-Hugoniot relation. We also 
plot the mean density of our Model la and lb as crosses, 
which are evaluated at the times depicted in Fig. [21 and 
plot the estimate of the critical mean density in equa- 
tion ([T5]) as a dotted line. If the mean density is larger 
than the critical density, a SN shock wave can gener- 
ate molecular cloud by piling up the shell. The results 
of simulations agree well with theoretical estimation of 
equation ([T6j) . From this figure, we expect that the for- 
mation of molecular clouds by piling up the two-phase 
shell is possible roughly when < 30 deg.. 

We conclude that the formation scenario of molecular 
clouds from the two-phase medium (HI clouds embed- 
ded in WNM) is expected to be a typical scenario. In 
this scenario, the fundamental building blocks of molec- 
ular clouds are not WNM but H I clouds consistent with 
the suggestions from observations of molecular clouds in 
Local Group galaxies (Fukui 2007; Blitz et al. 2007). 
Recently, we showed that a very strong shock wave that 
sweeps up two-phase gas cause turbulence and magnetic 
field amplification in the post shock medium (Inoue et 
al. 2009, see also Giacalone & Jokipii 2007). A similar 
phenomenon would also affect the physical conditions of 
molecular clouds formed from HI clouds as a result of 
shock sweeping that will be examined in a future paper. 

Note that the above speculation is obtained by assum- 
ing that the mean magnetic field strength in WNM is 
near equipartition to the thermal energy (B ~ 3 /iG). If 
the mean magnetic field strength is much larger than the 
equipartition, the direct formation of molecular clouds 
becomes easier, since the critical velocity derived in eq. 
(14) depends linearly on B. Also, the weak magnetic 
field would lead the direct formation of molecular clouds, 
if magnetic pressure is too weak to stop the isobaric con- 
traction of the shocked slab. 
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